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ABSTRACT 

We have mapped the central region of the HH 111 protostellar system in 1.33 
mm continuum, C 18 (J=2-l), 13 C0 (J=2-l), and SO (Nj = 5 6 - 4 5 ) emission 
at ~ 3" resolution with the Submillimeter Array. There are two sources, VLA 
1 (=IRAS 05491+0247) and VLA 2, with the VLA 1 source driving the HH 
111 jet. Thermal emission is seen in 1.33 mm continuum tracing the dust in 
the envelope and the putative disks around the sources. A flattened, torus-like 
envelope is seen in C 18 and 13 CO around the VLA 1 source surrounding the 
dust lane perpendicular to the jet axis, with an inner radius of ~ 400 AU (1"), an 
outer radius of ~ 3200 AU (8"), and a thickness of ~ 1000 AU (2"5). It seems to 
be infalling toward the center with conservation of specific angular m omentum 



rather than with a Keplerian rotation as assumed by I Yang et al.l (119971 ). An inner 
envelope is seen in SO, with a radius of ~ 500 AU (1''3). The inner part of this 
inner envelope, which is spatially coincident with the dust lane, seems to have a 
differential rotation and thus may have formed a rotationally supported disk. The 
outer part of this inner envelope, however, may have a rotation velocity decreasing 
toward the center and thus represent a region where an infalling envelope is in 
transition to a rotationally supported disk. A brief comparison with a collapsing 
model suggests that the flattened, torus-like envelope seen in C 18 and 13 CO 
could result from a collapse of a magnetized rotating toroid. 

Subject headings: stars: formation — ISM: individual: HH 111 — ISM: circum- 
stellar matter. 
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Introduction 



Stars are formed inside molecular cloud cores by means of gravitational collapse. The 
details of the process, however, are complicated by the presence of magnetic fields and angular 
momentum. As a result, in addition to infall (or collapse), rotation and outflow are also seen 
toward star-forming regions. In theory, a rotationally supported disk is expected to form 
inside a collapsing core around a protostar, from which part of the material is accreted by the 
protostar and part is ejected away. Observationally, however, when and how a rotationally 
supported disk is actually formed are still unclear, because of the lack of detailed kinematic 
studies inside the collapsing core. 



A rotationally supported disk has been seen with a r adius of 
(i.e., Class II or T Tauri) phase of star formation (see, e.g. 



500 AU in the late 



Simon et al.ll2000i ). Formation 



of such a disk must have started early in the Class phase. However, no such disk has 
been confirmed in the Class phase, probably because it is still too small (< 100 AU) for a 
detai led kinematic study with current instruments (see, e.g., Lee et al. 20061 ; Jprgensen et al. 



20071 ). In order to investigate how a rotationally supported disk is formed inside a collapsing 
core, we have mapped the protostellar system HH 111, which is in the Class I phase and 
is thus expected to have a bigger disk than that in the Class phase. In this system, a 
Class I source IRAS 05491+0247 is found deeply embedded in a compact molecular cloud 
core in the L1617 cloud of Orion, driving the HH 111 jet (IReipurth Il989h . The distance 



to the Orion Nebul a Cluster has been recently refined to be 400 pc ( jSandstrom et al.l 12007; 
Menten et al.ll2007l) . the same that is generally adopted for the Orion B cloud (see, e.g.. lGibb 
20081 ) with which L1617 is associated. The HH 111 system has traditionally been assumed 
to be at a distance of 450 pc, but w e here adopt this new distance, at which the source has 
a luminosity of ~ 20 (lReiDurthlE989h. A rot ating envelope has been seen around the 



source inside the core (jStapelfeldt fc Scovilld 119931 ) and may have infall motion toward the 
source ( Yang et al. 19971). The jet comp lex is highly collimated with a two-sided length of ~ 
6.7 pc (IReipurth. Bally, fc Devind 1 199 71 corrected for the new distance), a nd it produces a 
collimated molecular outflow around it ( Nagar et al. 1997 ; Lee et al. 200oh. The sourc e has 
also been detected with the VLA, associate d with a tiny radio jet (IReipurth et al.lll999l ). and 
a tiny dusty disk (with a radius of 30 AU) ([Rodriguez et al.ll2008l ) that can be the innermost 
part of a rotationally supported disk. Thus, this system is one of the best candidates to 
search for a rotationally supported disk and investigate how it is formed inside a collapsing 
core. Here, we present observations of the envelope in 1.33 mm continuum, 13 CO(J = 2 — 1), 
C 18 0(J = 2 — 1), and SO (Nj = 5q — 4 5 ) emission obtained with the Submillimeter Array 
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(SMA)0 JHo et al.ll200"3 ) 



2. Observations 

Observations toward the protostellar system HH 111 were carried out with the SMA 
on 2005 December 5 in the compact-north configuration using the 230 GHz band receivers. 
The receivers have two sidebands, lower and upper. With the local oscillator (LO) frequency 
set to 224.677 GHz, we observed the 12 CO(J = 2 — 1) line in the upper sideband, the 
13 CO(J = 2-1), C 18 0(J = 2-1), and SO (Nj = 5 6 - 4 5 ) lines in the lower sideband 
simultaneously. The 12 CO line traces mainly the outflow interaction and will be presented in 
a future publication. Combining the line-free portions of the two sidebands results in a total 
continuum bandwidth of ~ 3.9 GHz centered at the LO frequency (or A ~ 1.33 mm) for the 
continuum. Seven antennas were used, giving baselines with projected lengths ranging from 
8 to 53 m. The correlator was set up to give a velocity resolution of 0.28 km s -1 per channel 
for the 12 CO, 13 CO, and C 18 lines, and 0.56 km s -1 per channel for the SO line. Three 
pointings with a separation of ~ 28" were used to map the envelope and the outflow lobe 
extending to the west. Here, we only present the envelope region close to the source. 

The visibility data were calibrated with the MIR package, with Uranus as a flux cali- 
brator, quasar 3c454.3 as a passband calibrator, and quasar 0530+135 as a gain calibrator. 
The flux uncertainty is estimated to be ~ 15 %. The calibrated visibility data were imaged 
with the MIRIAD package. The dirty maps that were produced from the calibrated visibility 
data were CLEANed using the Steer clean method, producing the CLEAN component maps. 
The final maps were obtained by restoring the CLEAN component maps with a synthesized 
(Gaussian) beam fitted to the main lobe of the dirty beam. With natural weighting, the 
synthesized beam has a size of 3''5x3!'2 at a position angle (P.A.) of ~ —50°. The rms noise 
levels are ~ 0.22 Jy beam -1 in the 13 CO and C 18 channel maps, 0.15 Jy beam -1 in the SO 
channel maps, and 3.5 mJy beam -1 in the continuum map. The velocities of the channel 
maps are LSR. 



lr rhe Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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Results 



Our results are presented in comparison to a mosaic image based on the Hubble Space 
Telescope (HST) NICMO S ima ges ([Fell] 1.64 //m + H 2 at 2.12 /im + continuum) ob- 
tained by iReipurth et al.l (119991 ). which shows clear detection of two infrared sources, re- 
flection nebulae tracing the illuminated outflow cavity walls, and the jet in the system. 
The two infrared sources are named stars A and B, with star A being the IRAS source. 
They are also detected in 3.6 cm continuum by the VLA as the VLA 1 and 2 sources, 



respectively, at atpmn) = 05 h 51 m 46?25, S 



(2000) 



+02°48'29"5 and a 



(2000) 



05 h 51 m 46?07, 



<5(2ooo) — +02°48'30!'6 (IReipurth et al.lll999l ). These VLA positions, which are more accurate 



than the NICMOS positions, are adopted here as the source positions. In this region, the 
systemic velocity is assumed to be given by the line peak of the optically thin SO emission, 
which is 8.84±0.28 km s -1 LSR. Note that this valu e is higher than tha t deriv ed from the 
13 CO emission averaged over a much larger region by IReipurth fc Olbergi (Il99ll ). which was 
8.5 km s -1 . Throughout this paper, the velocity is relative to this systemic value. 



3.1. Continuum Emission 



A compact continuum emission is detected at 1.33 mm, with an emission peak well 
coincident with the VLA 1 position (Fig. []Ji). It has a total flux of 280±40 mJy, about 60% 
of th at measured for a much larger region in the single-dish observation (IReipurth et al. 
19931 ). The emission is spatially unresolved, mainly arising from around the VLA 1 source. 
Some faint emission is also seen extending towards the VLA 2 source. From the spectral 
energy distribution (SED) (Fig. [lb), it is clear that the emission at 1.33 mm is thermal 
emission from the dust in th e envelope and the pu tative disks around the VLA 1 and 2 
sources (IReipurth et al.l Il999l ; iRodrfguez et al.l 120081 ) . Assuming that the mass opacity is 
given by 



k v = 0.1 



cm 2 g 1 



10 12 Hz, 

the SED between A =3.06 mm and 25 /j,m can be fitted with (3 ~ 0.9 and a dust temperature 
of 41—64 K. With these fits, the emission is found to be optically thin at 1.33 mm. The total 
mass of the gas and dust associated with the envelope and the putative disks is estimated 
to be ~ 0.15-0.09 M . 
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3.2. Line Emission 

Unlike the continuum emission, C 18 0, 13 CO, and SO emission are all seen with a peak 
near the center of the dark ridge (i.e., the dust lane) seen in the HST image, which is ~ 0"5 
to the west of th e VLA 1 source (Fig. [2]). The C 18 emission shows similar structures to 



those seen in CS (jYang et al.lll997l ). It shows a north-south elongation roughly perpendicular 
to the jet axis, tracing a flattened envelope around the VLA 1 source. The C 18 emission 
also extends slightly to the west along the jet axis and southwest along the outflow cavity 
wall outlined by the reflection nebulae. The 13 CO emission is more extended, tracing both 
the envelope and the outflow cavity walls. The SO emission is more compact with a narrow 
waist around the source. It also extends along the outflow cavity walls in the southeast, 
southwest, and northwest directions. 

The dark ridge, where the line emissions all peak, ha s been suggested to be the birthplace 



of both the VLA 1 and 2 sources (IReipurth et al.lll999l ) and thus can be assumed to be the 



center of the flattened envelope. The C 18 0, 13 CO, and SO spectra toward the dark ridge all 
show a similar FWZI of ~ 6 km s _1 (Fig. [3]). The SO spectrum is rather symmetric with 
a single peak at the systemic velocity. The C 18 and 13 CO spectra show a double-peaked 
profile with a dip coincident with the SO peak, and with the blueshifted peak brighter than 
the redshifted peak (i.e., the blue asymmetry). Comparing to the single-dish observations 
toward a larger region, we find that there is missing flux in 13 CO and C 18 within ~ 1 
km s _1 of the systemic velocity as the large-scale flux is resolved out in our interferometric 
observation (Fig. H]). Note that, however, since missing flux tends to produce a dip with 
similar redshifted peak and blueshifted peak, an additional mechanism is still needed to 
produce the observed dips with the b lue asymme try. It could be self- absorption due to an 



inf all motion in the e nvelope (see, e.g.. lEvanslll999l ). Infall motion has already been suggested 



by lYang et al.l (119971 ) using CS observations. The 13 CO emission is optically thicker and thus 
shows a deeper and wider absorbing dip than the C 18 emission. The SO emission does not 
show a dip, likely because it is optically thin due to low SO abundance. 

In order to show the kinematics of the flattened envelope, the emission is divided into 
four velocity components: low redshifted and blueshifted emission with velocity within ~ 
1 km s _1 from the systemic, and high redshifted and blueshifted emission with velocity 
higher than 1 km s _1 from the systemic (Fig. [2]). At low velocity, the C 18 emission clearly 
shows a flattened structure with the blueshifted and redshifted emission extended ~ 7" (2800 
AU) to the north and south, respectively, from the dark ridge, as expected for a rotating 
envelope. The redshifted emission fits right in between the eastern and western reflection 
nebulae in the south. This is to be expected since the reflection nebulae trace the regions 
with less extinction and thus low density in the boundary of the envelope. Therefore, with 
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the boundary defined by those reflection nebulae, the envelope is found to have a thickness 
of ~ 2"5 (i.e., 1000 AU) and roughly constant with increasing distance from the source. Note 
that the blueshifted emission also traces the outflow emission in the west and the outflow 
cavity wall in the southwest. At high velocity, the emission shrinks toward to the source, 
arising from the inner part of the envelope. The blueshifted peak shifts slightly to the west 
and the redshifted peak slightly to the east. These shifts further suggest an infall motion 
in the envelope, with the blueshifted peak in the far end (west) of the envelope and the 
redshifted peak in the near end (east), as suggested in lYang et al.l (119971 ). These shifts are 
also partly due to an outflow component, with the emission extended to the west and along 
the cavity walls. In 13 CO, the blueshifted emission peaks to the northwest and redshifted 
emission in the southeast of the dark ridge even at low velocity. As seen in the spectrum, 
the lowest- velocity 13 CO emission, which traces the flattened envelope further out, is missing 
more than the C 18 emission, and thus revealing the inner part of the envelope, where the 
radial motion becomes prominent. In SO, the blueshifted and redshifted emission are also 
seen to the northwest and southwest of the dark ridge, similar to that seen in 13 CO. The SO 
emission is likely optically thin, probing directly the inner part of the envelope. 



4. A simple model for the envelope 

4.1. Velocity structure of the envelope 

The rotational velocity of the flattened envelope can be studied with the position- 
velocity (PV) diagrams centered at the dark ridge cut perpendicular to the jet axis. The 
13 CO and C 18 emission show similar kinematics (Figs. [5^, b) and can be studied together. 
The blueshifted emission is seen mainly to the north and redshifted emission mainly to the 
south, confirming a rotation in the envelope. At low velocity, the blueshifted and redshifted 
emissions also extend across the center to the south and north, respectively, suggesting that 
there is also a radial motion in the envelope. The rotational velocity is seen increasing 
toward the center and thus could be Keplerian or that conserving specific angular momen- 
tum. In the case of Keplerian rotation, the best fit to the velocity structures in these two 
line emissions together is = 1.3(R/ Rq)~ 1 ^ 2 km s" 1 with R = 400 AU (i.e., 1") , whic h 



is slightly faster than that found in CS at lower angular resolution by [Yang et al.l (119971 ). 
This Keplerian rotation corresponds to a stellar mass of ~ 0.8 M Q . Here, the plane of 
the flattened envelope is assumed to have an inclination of 10° away from the line of sight 



( iReipurth. Raga. fc Heathcotdll992l ). In the case of conservation of specific angular momen- 
tum, the best fit is = 1.6(RJ Ro)' 1 km s" 1 . With a reduced \ 2 — 0.07, it seems to fit 
better than the Keplerian rotation, which has a reduced \ 2 = 0.11 and decreases too fast 
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toward the center as compared with the observed. To confirm it, however, both observations 
at higher angular resolution and shorter uv spacings are needed. The latter is to retrieve 
the missing flux at low velocity. The SO emission traces only the inner part of the envelope. 
Its velocity structure is not well resolved and seems to have two components: a low-velocity 
component in the inner part with the velocity within ~ 1.5 km s _1 of the systemic, and a 
high- velocity component in the outer part with the velocity higher than ~ 1.5 km s _1 (Fig. 
[5b). The low- velocity component seems to have a radius of ~ 0''5 (200 AU) and a differen- 
tial rotation with the velocity increasing toward the center. The high-velocity component 
extends to R ~ 1.3" (~ 500 AU), connecting to that of the 13 CO emission. Its velocity, 
however, seems to increase with increasing distance from the center. It could arise from a 
solid-body (rigid) rotation with = 2.0(R/Rq) km s -1 . It could also arise from an outflow 
motion. 

The radial motion in the envelope can be studied with the PV diagrams cut along the 
jet axis (Figs. [5bl, e, f). The PV diagrams, however, are highly contaminated by the outflow 
emission, especially on the blueshifted side. Only part of the emission close to the center to 
within e.g., 2" (which corresponds to a deprojected distance of ~ 10" at an inclination of 
~ 10°) is believed to be from the envelope. The C 18 0, 13 CO, and SO emission within ~ 2" 
from the center seem to show similar velocity structure. The blueshifted emission is mainly 
in the west and the redshifted emission is mainly in the east. On the redshifted side, where 
the outflow contamination is less, the velocity seems to increase toward the center. This 
velocity structure is unlikely due to an outflow mot ion for which the velocity is expected to 



decrease toward the center (see, e.g. iLee et al.ll200ll ). It could arise from an infall motion in 



the envelope, as suggested . For a dynamica l collapse, the infall motion can be assumed to 
be a free-fall motion as in lYang et al.l (Il997l ). with the infall velocity v r = — v r o(R/ Ro) ~ l l 2 . 



Three values of v r0 , 2.4, 1.9, and 1.5, are compared to the C O, CO, and SO ve locity 



structures, with the second value corresponding to that obtained by lYang et al.l (119971 ). The 
corresponding stellar masses are 1.3, 0.8, 0.5 M G , respectively. The one with the least mass 
seems to give the best fit, matching the emission at the lowest velocity in the redshifted 
part where the contamination of the outflow emission is less. However, observation at higher 
angular resolution is really needed to confirm it. 



4.2. Density structure of the envelope 



In order to obtain roughly the de nsity distri b ution of the flattened envelope, a simple 
disklike envelope model as described in lLee et al.l (120061 ) is compared to the C 18 emission, 
which suffers less missing flux than the 13 CO emission. In this model, the envelope has a 
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free-fall motion and a rotation with conservation of specific angular momentum, with the 
best-fit parameters given above. It has an inner radius of i?; n , an outer radius of R ut, 
and a constant thickness of H, in cylindrical coordinates. The number density of molecular 
hydrogen is assumed to be given by 



n = n (— ) 



(2) 



where no is the density at Ro and p i s a power- l aw index assume d to be —1.5, as in many 
theoretical infalling models (see, e.g., IShu! 119771 ; iNakamural l2000l ) . The temperature of the 
envelope is uncertain and assumed to be given by 



(3) 



where To is the temperature at Ro and q is a power-law index. The abundance of C 18 Q relative 
to m olecular hydrogen is assumed to be constant and given by 1.7 x 10~ 7 ( IFrerking et al. 
19821 ). In the model calculations, radiative transfer is used to calculate the emission, with 
an assumption of local thermal equilibrium. For simplicity, the line width is assumed to be 
given by the thermal line width only. The channel maps of the emission derived from the 
model are convolved with the observed beam and velocity resolution, and then used to make 
the integrated map, spectrum, and PV diagram. 

There are six parameters in this model: R m , R on t, H, no, To, and q. In our fits, we have 
R in ~ 400 AU (~ 1"), .Rout ~ 3200 AU (~ 8"), H ~ 1000 AU (i.e., 2%), n ~9x 10 6 cm" 3 , 
Tq ~ 40 K, and q ~ —0.6. Here T is assumed to be the lower end of the dust temperature. 
As can be seen in Fig. 0, this model can roughly reproduce the envelope structure, spectrum, 
and PV diagram of the C 18 emission, except for the emission at high-redshifted velocity. 
The position of the two peaks in the line profile can also be roughly matched with the 
assumed infall velocity. Note that, however, the temperature power-law index is required to 
be steeper here in order to reproduce the observed dip at the systemic velocity, as compared 
to those (q y —0.4 ) obtained by fitting other envelopes in single-dish observations (see, e.g., 
Hogerheijdd l200ll ). It could be because the dip is partly due to the missing flux of the 
large-scale gas and the absorption by this large-scale gas as well. 



The mass of the flattened envelope, which is given by 

" R out ft 



2mnnoH 



R\) 



1.5 



2nRdR 



(4) 



is ~ 0.2 M Q , consistent with that found by lStapelfeldt fc Scovilld (119931 ). This mass is higher 
than that derived from the continuum emission, which is more compact. In our simple model, 
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the infall rate, which is given by 

M{R) = 2irRHn\v R \ • 2m H 

~ 4 x 10~ 5 A" 1 M yr- 1 (5) 
ito 

has a mean value of ~ lO~ 5 M yr _1 , averaged over the radius. The accretion rate, if assumed 
to be the mean infall rate, would result in an accretion luminosity of L acc = GM*M acc /R* ~ 
37 L Q , assuming a stellar mass of M* ~ 0.5 M Q (as derived from the infall velocity) and a 
stellar radius of R* ~ AR ^. This luminos ity, however, is about twice the luminosity of the 



source, which is ~ 20 L Q (jReipurthlll989l ). In addition, the accretion time would be only ~ 
5xl0 4 yr for forming a source with a mass of 0.5 M Q , shorter than that for a typical Class 
I source. Thus, the accretion rate could be smaller than the infall rate. 



4.3. Summary - A Collapsing torus 

The envelope seen in C 18 and 13 CO can be considered as a thick torus around the 
narrow dust lane (i.e., the dark ridge), with an inner radius of 400 AU (1"), an outer radius 
of 3200 AU (8"), and a thickness of 1000 AU (2'/5). It has a differential rotation that seems 
to be better fitted by a rotation w ith conservation o f specific angular momentum than by a 



Keplerian rotation as assumed by lYang et al.l (119971 ). The envelope likely also has an infall 



motion, with the blueshifted emission in the far end and the redshifted emission in the near 
end, showing a double-peaked line profile with a blue asymmetry and an absorption dip at 
the systemic velocity. Therefore, the envelope seems to be infalling (contracting) toward the 



(Ohashi et al 


1997 




Lee et al. 


2005) 



phase (e.g., iLee et al.ll2005l ). With the derived rotation and infall velocity structures, the 
rotation velocity is equal to the infall velocity at R ~ 1"5 (600 AU), suggesting that the 
centrifugal force starts to dominate in the inner edge of the envelope. 



5. Discussion 



5.1. Magnetic braking and SO disk? 



The inner envelope is seen in SO as in HH 212 (ILee et al.ll2006l ). It has an outer radius 
similar to where the centrifugal force starts to dominate. Therefore, a rotationally supported 
disk is expected to form here. Such a disk may indeed have formed in the inner part of this 
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inner envelope, which is seen with a differential rotation, spatially coin cident with the narrow 
dark ridge that could be interpreted as a dense circumstellar disk (IReipurth et al.l Il999l ) . 
The outer part of this inner envelope, however, may have a rotation velocity decreasing 
toward the center. If this is the case, the angular momentum might have been efficiently 
removed from there. This can be achieved with the presence of a magnetic field. It has been 
found that a magnetic field can enforce solid-body rotation because magnetic braking can 
eliminate differential rotation on a timescale much shorter than even the free-fall timescale 
rtMouschovias & Cioled[l~9~9~9[ ). 



5.2. Comparing with a Theoretical Model 



Currently, the most popular model of star formation is Shu's model (derived from IShu 



19771 ) . One recent ve rsion of this model assumes the collapse of a magnetized singular isother- 
mal rotating toroid (jAllen. Li. fc Shull2003l ; iMellon fc Lil 120081 ). In this model, a pseudodisk 
(a nonrotationally supported thick torus) is formed around the source. In the inner edge 
of the pseudodisk, the rotation speed drops steadily with decreasing radius, as a result of 
efficient braking by the strong magnetic field. 

In our observations, it is possible that the flattened envelope seen in C 18 and 13 CO is 
a pseudodisk resulting from the collapse of a magnetized rotating toroid toward the center. 
It is also possible that the outer part of the SO envelope is the inner edge of the pseu- 
dodisk where the rotation velocity decreases with decreasing radius due to efficient magnetic 
braking. However, further observations are really needed to confirm the rotation velocity 
structure in the flattened envelope and to tell whether the structure of the infalling mate- 
rial is magnetically-guided. In a ddition, a slow outflow , which is expected when there is a 
twisting of magnetic field lines (jAllen. Li. fc Shu! 120031 ). has not yet been identified in our 
observations. 



In their model, a rotationally support ed disk is not able to form around the sou rce due 
to the efficient loss of angular momentum (jAllen. Li. fc Shull2003t IMellon &: Lill2008l ). In the 
HH 111 system, such a disk must have already formed around the source to launch the jet. 
In addition, the innermost pa rt of the disk has been seen as a tiny du sty disk (w i th a r adius 
of 30 AU) around the source (jRodriguez et al.ll2008l ). As suggested by lShu et al.l (120061 ) . one 
possibility is to have an ohmic resistivity in the model, so that the magnetic field can be 
dissipated efficiently around the source. 
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5.3. Gas depletion? 



The continuum emission peaks at the VLA 1 position while the line emissions all peak 
around the dark ridge, suggesting that the dust and gas peak at different positions. It could 
be because there is a lack of gas toward the continuum peak. It has been suggested by 
Reipurth et al.l (119991 ) that the VLA 1 and VLA 2 sources were originally formed in the dark 
ridge inside the envelope but later ejected from there. Thus, the continuu m emission peak 



could trace the inner disk carried away by the VLA 1 source, as suggested in lReipurth et al. 



(119991 ) . If this is the case, there could be a lack of gas (depletion) toward the inner disk. 



6. Conclusion 

We have mapped the central region of the protostellar system HH 111 in 1.33 mm 
continuum, C 18 (J=2-l), 13 CO (J=2-l), and SO (Nj = 5 6 -4 5 ) emission at ~ 3" resolution 
with the Submillimeter Array. There are two sources, VLA 1 (=IRAS 05491+0247) and 
VLA 2, with the VLA 1 source driving the HH 111 jet. Thermal emission is seen in 1.33 mm 
continuum tracing the dust in the envelope and the putative disks around the sources. A 
flattened, torus-like envelope is seen in C 18 and 13 CO around the VLA 1 source surrounding 
the dust lane perpendicular to the jet axis, with an inner radius of ~ 400 AU (1"), an outer 
radius of ~ 3200 AU (8"), and a thickness of ~ 1000 AU (2"5). It seems to be infalling toward 
the center with conserv ation of specific a ngular momentum rather than with a Keplerian 



rotation as assumed by lYang et al.l (119971 ). An inner envelope is seen in SO, with a radius 
of ~ 500 AU (1''3). The inner part of this inner envelope, which is spatially coincident with 
the dust lane, seems to have a differential rotation and thus may have formed a rotationally 
supported disk. The outer part of this inner envelope, however, may have a rotation velocity 
decreasing toward the center and thus represent a region where an infalling envelope is 
in transition to a rotationally supported disk. A brief comparison with a collapsing model 
suggests that the flattened, torus-like envelope seen in C 18 and 13 CO could be a pseudodisk 
(a nonrotationally supported torus) resulting from a collapse of a magnetized rotating toroid. 



We thank the SMA staff for their efforts in running and maintaining the array. We 
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the NASA Astrobiology Institute under Cooperative Agreement No. NNA04CC08A issued 
through the Office of Space Science. 
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Fig. 1.— (q) The continuum emission overplotted on the HST NICMOS image from 
Reipurth et al.l (119991 ) . The asterisk marks the position of the IRAS source (= star A = 
VLA 1) and the cross marks the position of star B ( = VLA 2). The contours go from 
4 to 74 a with a step of 10 a, where o = 3.5 mJy beam -1 . The blue and red arrows in- 
dicate the orientations of the blueshifted and redshifted parts of the jet, respectively, (b) 
SED tow ard the central r e gion. The filled squar e is from our obse rvation. The diamonds 
are from iRodriguez et al.l (120081) . the star from lYang et al. Jl997h. t he op en circles from 



Stapelfeldt fc Scoville! Jl993f ). the filled triangles from iReipurth et ali (|l993h. the op en tri 



angles from the IRAS catalogue, and the open squares from lReipurth &: Olbergi (Il99ll ). The 
dashed and solid lines indicate the fits with a dust temperature of 41 K and 65 K, respectively. 
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Fig. 2.— C 18 (first column), 13 C0 (second column), and SO (third column) emission 
contours on the HST NICMOS image. The asterisk marks the position of the IRAS source 
(= star A = VLA 1) and the cross marks the position of star B ( = VLA 2). "x" marks 
the central position of the dark ridge and the dashed line indicates the equatorial plane 
perpendicular to the jet axis. (Top row) Total emission integrated from ~ —3 to 3 km s -1 . 
The contours go from 3 to 18 a with a step of 3 a for C 18 0, from 5 to 47 a with a step of 
7 a for 13 CO, and from 3 to 24 a with a step of 4 o for SO, where a = 0.3 Jy beam -1 km 
s -1 . (Middle row) Low blueshifted and redshifted emission integrated from ~ — 1 to km 
s -1 and from to 1 km s -1 , respectively. The contours go from 3 to 12 a with a step of 3 a 
for C 18 0, from 5 to 29 a with a step of 6 a for 13 CO, and from 3 to 21 o with a step of 6 a 
for SO, where a = 0.11 Jy beam -1 km s -1 . (Bottom row) High blueshifted and redshifted 
emission integrated from ~ —3 to —1 km s -1 and from 1 to 3 km s -1 , respectively. The 
contours go from 3 to 12 a with a step of 3 a for C 18 0, from 5 to 40 a with a step of 7 a for 
13 CO, and from 2 to 8 a with a step of 2 a for SO, where a = 0.16 Jy beam -1 km s -1 . 



-16- 




-3 
Velocity offset from systemic (km/s) 

Fig. 3. — C 18 (thin black), 13 C0 (thick black), and SO (thick gray) spectra toward the 
dark ridge seen in the NICMOS image. The 13 CO spectrum has been multiplied by a factor 
of 0.3. 
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Comparison of our C 18 and 13 CO spectra (in gray) with those obtained with the 
30 m telescope (in black) (ICernicharo fc Reipurthlll996l ). 
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Fig. 5. — Position-velocity (PV) diagrams centered at the dark ridge. (Top row) PV cuts 
perpendicular to the jet axis. In (a) and (b), the dots with error bars represent the peak 
positions determined at each velocity channel. Solid curves are derived from the rotation with 
conservation of specific angular momentum. Dashed curves are derived from the Keplerian 
rotation. In (c), the solid lines are from the linear fit to the outer part of the SO envelope. 
(Bottom row) PV cuts along to the jet axis. Solid, dashed, and dot-dashed curves are derived 
assuming a free-fall motion v r = —v^R/Rq)' 1 ^ 2 , with v r o = 2.4, 1.9, and 1.5, respectively. 
The contours start at 2 a and have a step of 2 a, where a = 0.2 Jy beam" 1 for C 18 and 
SO and o = 0.4 Jy beam" 1 for 13 CO. 
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Fig. 6. — A simple model fitting to the C 18 emission, spectrum, and PV diagram cut 
perpendicular to the jet axis. See the text for model description. Black contours and 
spectrum are from model. 



